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SUMMARY 


An investigation is made to obtain some information concerning the 
nature and magnitude of the errors involved in computing the landing 
stresses for large and relatively flexible aircraft when several simpli- 
fying assumptions are made. An exact solution is made for the landing 
stresses of a simplified structure and is compared with several approxi- 
mate solutions made when the simplifying assumptions are used. 

The simplified structure investigated consisted of a uniform beam 
for the wing, a concentrated mass for the fuselage, and an undamped 
linear spring for the landing gear. This structure was considered to be 
in uniform translation until the landing gear touched the ground. The 
subsequent motion was computed by using operational calculus in conjunc- 
tion with standard beam theory. In general. It was found that, for 
moderately flexible landing gears, the neglect of the effect of structural 
elasticity in computations of strut forces or of the acceleration of the point 
of attachment of the landing gear and then the computations of stresses by 
considering the structure to be elastic led to small conservative errors } 
whereas the neglect of structural elasticity in computing wing stresses 
from the strut force or acceleration of the point of attachment of the 
landing strut led to unconservative errors of appreciable magnitude. 

This result suggests that a satisfactory treatment of the landing problem 
may possibly be obtained from an analysis which assumes that in landing 
the aircraft Is an elastic structure subject to the forces or accelera- 
tions found in a drop test In which a rigid mass Is used. 


INTRODUCTION 


When an aircraft lands, the vertical component of its velocity Is 
rather suddenly reduced to zero. This sudden change in motion of the 
aircraft gives rise, to stresses within the structure which may become 
large and even destructive as the size and weight of airplanes increase 
and the design load factor decreases. The design of large airplanes 
should therefore consider the effect of a severe landing on the wings, 
fuselage, tail surface, landing-gear struts, and other elastic parts of 
the airplane structure. 
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When the stiffness of an airplane structure is relatively large 
compared with the stiffness of the landing gear, all parts of the 
airplane are subject to essentially the same acceleration during landing 
and the stresses in all parts of the airplane are therefore readily 
computed by methods of statics. When the stiffnesses of the landing 
gear and the airplane structure are comparable in magnitude, the various 
parts of the airplane have different accelerations and the problem 
becomes much more involved. The calculation of the stresses is a 
complex problem even when the equations involved are purely linear j that 
is, when the internal forces are proportional to the deformations. Xt 
is much more complicated in the actual aircraft because of the nonlinear 
characteristics of the landing gear. A number of approximate methods of 
computing landing stresses have consequently been used. 

Because of the practical importance of the landing problem, it is ofi 
some interest to determine the nature and magnitude of the error involved 
in various approximate-methods that- have been used. For this purpose, 
an exact solution is made to determine the landing stresses for a highly 
simplified structure in which a uniform beam is used to represent- the 
wing, a rigid mass the fuselage, and a simple spring the landing gear. 

The stresses in the wings excited by the landing impact are computed 
by operational calculus in conjunction with the standard engineering 
vibration theory of beams. The results are compared with the results 
found by a number of approximate methods. The analytical treatment of 
the exact and approximate solutions are given in appendixes. 


SYMBOLS 


E modulus of. elasticity 

y density of wing material in units of weight 

c velocity of sound in wing material 

g acceleration of gravity 

L semi span of wing 

I moment of inertia of cross section of wing about 

neutral axis 

A cross-sectional area of_wing 

p radius of gyration of cross section of wing 

x coordinate along wing measured from root 

y distance from neutral axis of wing to - any fiber 
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t 


m mass of wing .(seijiispan) 

M concentrated, mass (one-half of fuselage mass) 

S spring stiffness 


t 

P 

n 


'n 


time, zero at beginning of impact 
operator 

integers 1 , 2 , 3 , and so forth designating a 
particular mode of vibration 

nth positive root of transcendental equation 
associated with a given type of vibration 



“n 


undamped natural angular frequency of nth mode, radians 
per second 


v 

®c 

COQ 

CD 

S 

w(x, t) 

a(x, t) 
ct(x, y, t) 

t(x, t) 


vertical velocity of aircraft prior to impact 
natural frequency of fundamental mode of a cantilever. 


radians per second 




natural frequency of fuse lags -spring system, radians 
per second ^ || ^ 

natural frequency of airplane with wing rigid, radians 

per second ( \ — § — | 

IV M + m j 

deflection, relative to root position at t “ 0 , 
of wing at station x and time t 

acceleration of wing at station x and time t 

bending stress in wing at station x, distance from 
neutral axis y, and time t 

average shear stress over cross section of beam at 
station x and time t 


Ajj bending -stress coefficient 

A maximum bending-stress coefficient obtained from 

first three modes with proper regard to phase 
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B n shear-stress coefficient 

B maximum shear-stress coefficient obtained from first - 

three modes with proper regard to phase 


RESULTS AND DISCUSSION 


Exact solution. - In order to obtain pertinent information on the 
problem of landing impacts, an exact solution is made of the landing 
stresses of a highly simplified structure. In this simplified structure 
(fig. l) , a uniform beam of mass m was used to represent the wings of 
the airplane, a rigid mass of magnitude M to represent the fuselage, 
and a simple spring ot stiffness S was substituted for the landing 
gear. The exact analytical treatment giving the equations for frequencies, 
deflections, accelerations, strut force, bending stresses, and shear 
stresses is presented in appendix A. The maximum root bending stress 
that results from impact (gravity not included) is shown in this appendix 
to be given by the equation 

cr = A - Je 
c p 

where A is a dimensionless coefficient dependent on the physical 
parameters of the structure. As can be Been, the maximum stress is 
directly proportional to the velocity of descent; 

In figure 2 the coefficient A is given for several values of the 
ratio of fuselage mass to wing mass as a function of the ratio 

In this ratio, co c is the fundamental frequency of the wing as a canti- 
lever and cd b is the frequency of the airplane when the wing is 
considered rigid. Low values of the frequency ratio correspond to a 
flexible landing gear and the corresponding induced stresses are rela- 
tively small but become larger as the landing gear becomes stiff er. The 

limiting case of a rigid landing gear m was investigated in 

reference 1, damping being taken Into account. The results showed that- 
damping eliminates the higher frequencies much faster than the lower ones 
so that only the lower modes might be expected to contribute to the 
maximum root bending stress. In the computation of the curves shown in 
figure 2, only the first three modes were considered, with proper regard 
being given to phase. (See appendix A for the stress that is associated 
with each mode.) On this basis A Is approximately equal to 2.8 in the 
case of a rigid landing gear. 
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Appro xima te solutions. - The exact solution Just discussed vas obtained 
■by solving the equations of motion directly. In the approximate solutions, 
the problem is broken arbitrarily into two parts or stages as follows: 

Stage 1: Determination of the strut reaction or of tfie acceleration 

of the points of attachment of the landing gear. 

Stage 2: Computation of the stresses in the airplane structure by 

use of one of the quantities obtained in stage 1. 

In both stages of this approach the structural elasticity muBt be 
properly taken into account if the correct solution is to be obtained. 

In the approximate solutions the effects of the structural elasticity 
are neglected in one or both stages. Five approximate solutions are 
given in appendix B. For convenience in discussion, the approximate 
methods are identified herein as follows: 

Method A - Structural elasticity neglected in both stage 1 and stage 2. 

Method B - Structural elasticity considered in determining reaction 
(stage 1) but neglected in stage 2. 

Method C - Structural elasticity considered in determining accelerations 
(stage 1) but neglected in stage 2. 

Method D - Structural elasticity neglected in determining reaction 
(stage 1) but considered in stage 2. 

Method E - Structural elasticity neglected In determining accelerations 
(stage 1) but considered In stage 2. 

Method F - Statistical approach of Biot and Bisplinghoff (reference 2) j 
structural elasticity considered In Btage 2. 

The simplest calculation, of course, results from use of Method A, 
which neglects the structural elasticity altogether. The acceleration 
of all parts of the structure is then assumed to be equal to the acceler- 
ation measured In a drop test in which a rigid mass equal to the mass of 
the airplane without the landing gear is used. The stresses are obtained 
by statics, from a wing loading obtained by multiplying the mass distri- 
bution by the acceleration found In the drop test. The coefficients for 
maximum root bending stress obtained by method A are shown in figure 3 . 

This curve Is Independent of the mass ratio M/m. For comparison the 

M 

exact solution for — = 2 Is also shown. This mass ratio Is used for 

m 

all of the succeeding comparisons . 

The physical assumptions of the two approximate methods in which 
the structural elasticity Is taken Into account in stage 1 but ignored 
In stage 2 are shown schematically in figure 4. The results found by 
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method B , where the true strut reaction (strut force given "by exact 
solution) is applied to a rigid airframe, and by method C, where the 
rigid airframe is subject to the true acceleration (acceleration given 
by exact solution) at the points of support are shown in figure 3* 


It appears from figure 3 that the three approximate methods, which 
neglect the effects of elasticity in the second stage, are unconservative 

0)_ / cn s \ 

for — 2<1 markedly so near — = 0.5 but approach the correct results 

“>c V ®o / 

as the ratio approaches zero. For aircraft with extremely stiff landing 


gears » 1 J, the approximate methods are highly conservative $ they 

predict infinite stresses when a rigid landing gear is used. The fact 
that curves obtained by methods B and C, which represent approximate 
methods in which structural elasticity is taken into acoount in the first 
stage but neglected in the second, agree much better with the curve 
obtained by method A, which neglects structural elasticity altogether, 
than with the exact solution suggests that the neglect of structural 
elasticity in stage 2 is much more serious than in stage 1. 


Two methods are then resorted to in whioh the effects of structural 
elasticity are neglected in stage 1 of the analysis but are properly 
taken into account in stage 2 (methods D and E) . The physical assumptions 
made are indicated schematically in figure 5 • In these methods the strut 
reaction and the acceleration are determined in a drop test in which a 
rigid mass is used. The resulting stresses are then computed with due 
regard for the elastic response of the structure. The results found are 
compared with the exact solution in figure 6. The curve shown for 
method D is for the strut -react! on method and the. curve for method E 
is for the acceleration method. The curves are cut off when o> B /c u c is 

about 0.7, since for higher ratios the force or accelerations obtained 
from a drop test with the simple undamped spring will give rise to 
resonance effects having very little relation to the actual landing 
problem. It appears from figure 6 that the two methods whioh neglect 
structural elasticity only in stage 1 (that is when determining strut- 
forces or the accelerations of the points of support) are conservative and 
are subject“to only small errors. 

A somewhat different method of handling the landing problem is the 
statistical approach developed by Biot and Bisplinghoff in reference 2. 

In this method (method F),the time history of-the landing impact is 
assumed to be independent of the elastic- properties of the structure, so 
that it may perhaps be classed with the methods which neglect structural 
elasticity in stage 1. A number of other approximations are also involved 
for the sake of simplifying the analysis and including a wide variety of 
landing conditions. Among these approximations are: 
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(a) The impact force is characterized "by only two parameters, the 
maximum value and the duration' of the force; thus, the detailed history 
of the collision is not considered. 

(h) The maximum stress in the first mode is obtained by assuming, 
in effect, that of a whole class of typical force-time histories having 
prescribed values for these two parameters; the force-time history that 
applies is the one which leads to the highest stress . 

(c) Similar assumptions are made in obtaining stresses in higher 
modes. This procedure results, effectively, in the assumption that, of 
the force-time histories used to determine the envelope response curve, 
a different one may apply to each mode . 

(d) The maximum stress is found by adding together the maximum 
stresses found for the first three modes without regard to phase. 

All of the approximations discussed in connection with the Biot and 
Bisplinghoff method are conservative except the use of only three modes. 
The restriction to three modes, which characterizes also the curves for 
the exact solution of the present paper, would be unoonservative in the 
undamped case; the airplane, however, is subject to a large amount of 
damping and no real unconservatism is likely to result. The expectation 
of a conservative result for the Biot and BiBplinghoff method (method F) 
is verified in figure 6. This method, which makes use of a number of 
conservative simplifying assumptions, appears in some cases to overestimate 
the stresses by a factor of almost 2. 


CONCLUDING REMARKS 


The problem of computing the landing stresses for a large and rela- 
tively flexible aircraft is so complex that most investigations are based 
on simplifying assumptions. The present paper constitutes an attempt to 
obtain some information concerning the nature and magnitude of the errors 
in these assumptions by solving the landing problem exactly for a simpli- 
fied structure and comparing the results with solutions to the same 
problem obtained by use of the simplifying assumptions. 

The simplified structure investigated consisted of a uniform beam for 
the wing, a concentrated mass for the fuselage, and an undamped linear 
spring for the landing gear. This structure was considered to be in 
uniform translation until the landing gear touched the ground. The 
subsequent motion was computed by using operational calculus in conjunction 
with standard beam theory. 

In most of the approximate treatments that have been proposed the 
problem is arbitrarily broken into two parts in the first of which the 
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strut reaction or the acceleration of the point— of attachment to the 
landing gear is determined (stage 1), and in the second of which the 
stresses resulting from the applied force or acceleration are calculated 
(stage 2) . The various approximate methods investigated either neglect 
altogether or treat only approximately the effects of structural 
elasticity in one or both stages. 

In general, for moderately flexible landing gears, the neglect- of 
the effects of structural elasticity in stage 2 was found to be more 
serious them the corresponding neglect- in stage 1. Such neglect in 
stage 2 led to unconservative errors of appreciable magnitude j in stage 1 
it led to errors which were smaller and on the conservative side . A 
statistical approach proposed by Biot~and Bisplinghoff -was found to be 
always conservative and to have in some cases a safety factor nearly equal 
to 2. 


The conclusions Just stated were based primarily on the analysis of ' 
the behavior of the simplified structure studied for a ratio of fuselage 
mass to wing mass equal to 2. The results, however, are essentially the 
same when the ratio is 1/2 or 5* It therefore appears reasonable to 
expect that the conclusions of-this paper have general validity as 
applied to conventional aircraft. In addition, the results suggest that 
a satisfactory treatment of the landing problem may possibly be obtained 
from an analysis which assumes that, in landing, the aircraft is an 
elastic structure subject to the forces or accelerations found in a drop 
test in which a rigid mass Is used. 


Langley Memorial Aeronautical Laboratory 

National Advl sory Committee for Aeronautics 
Langley Field, Va., January 30, 1948 



NACA TN No. 1584 


9 


APPENDIX A 
EXACT SOLUTION 


General analysis .- In order to make the problem of computing the 
landing stresses of an airplane susceptible to accurate solution without 
an inordinate amount of labor it is necessary to idealize the structure. 

The simplified structure used in the present analysis to represent the 
airplane is shown in figure 1. The airplane is; considered to be falling 
with a constant velocity v until the bottom of the spring is suddenly 
stopped by contact with the ground. This disturbance gives rise to 
oscillations in the beam governed by the differential equation (reference 1) 


Ep 2 




= 0 


(Al) 


Previous analyses (references 1 , 3 > ancL 4) have treated special cases of 
oscillations of a cantilever team due to impact, with an internal damping 
term included. Experience indicates, however, that damping haB only a 
slight effect upon the terms that are significant and therefore damping 
is neglected in the present analysis. 


With the use of the notation c 2 = ^ and the operational 

notation p = equation (Al) may be written as the ordinary fourth- 

ot 

order differential equation 




(A2) 


The general solution of this equation is 


where 


w = P cosh 9 ~ + Q sinh 0 — + B sin 0 — + S cos 0 £ (A 3 ) 

lx L L 



The coefficients P, Q,, E, and S are to be evaluated from the 
boundary conditions, which in this case are 
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0 



The last - boundary condition expresses the equilibrium of the forces that 
act on the mass M. The expression in "brackets on the right hand is the 
change in length of the spring and, when multiplied "by S, gives the force 

exerted "by the spring on the mass. The term /*(v - v/) dt indicates 


that the motion of- uniform translation at- the "bottom of the spring is 
suddenly stopped at t = 0, and the term w ( x _o) is the displacement of 

the root and is equal to the displacement of the "bottom of the spring 
for t < 0. The oscillations set up when the bottom of the spring is 
suddenly arrested from uniform translation would "be the same as if the 
bottom of the spring were suddenly set in uniform motion with the system 
initially at rest. The uniform- velocity term may therefore be omitted 

and the last boundary condition becomes, if v£) dt is replaced by 

v/ 

the operational form - 

P 



With the application of the boundary conditions to equation (A3) there » 
is obtained a set of- four nonhomogenous equations in terms of the four 
coefficients P, Q, P, and S. These equations are solved for the 
four coefficients and the equation for velocity may then be written. 

The operational form for the velocity (that induced in the beam when 
the bottom of the spring is suddenly set in motion) is found to be 


pw = 




(a4) 
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where 


3r(e^ = (l + cos 0 cosh 0)^cosh 0^ + cos 0 ^ 

+ sin 9 sinh Swoosh 021 _ C os 0^ 

- (cosh 0 sin 0 + sinh 0 cos 0) (sirih 0— - sin 0— j 

\ L L/ 


z -S. 


av 


CDu 

1 + l(l + cos 0 cosh 0 ) + ^(cosh 0 sin 0 + sinh 0 cos 0) 

p ^ 


and 


m 0=\/| 


Interpretation of equation (a 4) "by the Heaviside expansion theorem and 
addition of the constant velocity v gives for the total velocity 


MX’J'l = V - v/ 




+ V 


V 


“o 2 


cos 


“a* l 


(A5) 


where 0 n is the nth positive root of the equation 

2 N 


OV 


“n 


+ 1 cosh sln 9 n + slnh 0 n cos e n = Q 


M 


© n (l + cos 0 n cosh 0 n ) 


(A 6) 


and 


“n “ PC 72~ 


®n 


F ( 9 ^)‘ F ( e r) 


(with 0 replaced "by 0 n ) 
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A = [ 1 - ] -2 (sinh An cos ^ - cosh 6^ sin S^) 

“n 

+ 2^1 + cos 0 n cosh 9 ^ + 

+ - — ( cosh 0 sin 9^ + sinh 9 cos 6 \ 

2 0 n V n n n nj 

f 2 

co- 


rn 

Si 


cos 9 cosh 9 
n n 


2^1 + cos 0 n cosh 0^ 


m / . * 

- [ cos 9„ + cosh 0 A 
M \ n n / 


wq r 

1 + 3 (l + cos 0 n cosh 0 n V 

“n* ' 


03u *3 

0 qT tti "L 

The term — - may he transformed into the form =a~ - — r and thus 

o^ 2 EI M 9 n h 

equation (a 6) , which defines the root 0 n , may he written 


1 - 


Sir m 1 


EI M a 4 1 M 
°n 


_ cosh 0 sin 0 „ + sinh 0 „ cos 0 ,, 
, m n n n n 


= 0 


0 n (l + cos 0 n cosh 0 n ) 


Solution of this equation for SL^/EI gives 

SL 3 _ 0 4 M . 0 n 3 ( c °e h 0 n sln 0 n + sinh 9 n coa d n) 

^ ^ ■■■!■■ — — — — — — ■ ■- ■ - m A m 

TPT m 


EI 


1 + cos 0 cosh 0„ 
n n 


(A7) 


This equation (or any of its previous forms) is the characteristic 
frequency equation of the team-mass -spring system. A graphical represen- 

OT 3 

tation of this equation Is shown In figure 7. in which — — is plotted 
against 0 for values of S=0, 2, 5, 10, and 50. The values of 0 
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O 

corresponding to given values of* SL /EH and M/m are the roots 0 n 
which, characterize the modes of vibration. 


Integration of equation (A 5 ) with respect to time with the 
condition (w)^_^ = 0 gives for the. deflection 


w(x. 


CO 



n=l 



sin o^t £ 


(A8) 


From equation (A 5 ) for velocity and equation (a 8) for deflection the 
complete behavior of the idealized structure after landing may be found. 
The quantities of chief interest are the maximum bending stress, the 
maximum shear stress, the accelerations, and the force in the spring. 


Maximum bending stresses .- The maximum bending stresses cr(x, y, t) , 
at any fiber distance y from the neutral axis, occur at the root and 
are given by the equation 


ct(x, y, t) = Ey 


*2 
o w 


x=0 


= E 


v 7 


5 P Z_ A n sin o t / 
K n=l u 


(A 9 ) 


where 


ox. 


A =4 — 
■^•n ^ 2 


“0 2 


sin 0 n sinh 8 n + cos 0 n cosh 0^ 


1 + 3 — j(l + cos 0 n cosh 0 n ) + S (cos 0 n + cosh 0 n ) £ 

“n ' ‘ V 


Equation (A 9 ) may be written in the form 


cr * E c “ ( A 1 sil1 + ^2 sin + A 3 sin “3^ + • • • ) (AlO) 
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In the expression In parentheses each term may he thought-of as the 
contribution of a particular mode to the bending stress. The expression 
In parentheses may be denoted by A and equation (A10) becomes simply 


a 


= A - - E 
c p 


The values of A^, A^, and are plotted against 

for 2 , and 5 in figure 8 , and the maximum value of A found to 

m 2 

occur in the initial cycles of... vibration by use of the first- three terms 

of the series is plotted against 03 B y5n c for — = i, 2, and In figure 2 . 
The q uan tity co B is the frequency of the system If the ving were rigid: 


03 


s 



0 ^ 



m 

M 


The quantity 03 c is the frequency of the wing if the spring were 
infinitely stiff: 


a> c = pc 


3-52 

l2 


This equation represents the fundamental cantilever frequency. (See 
reference 1.) The ratio <jo b /cc c Is related to the parameters used in 

frequency equation (equation (A7)) by the relation 



1 

12.4(l + S) 



Maximum shear stress .- The maximum average shear stress t(x, t) 
occurs at the root and Is given by the equation 
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In figure 9 "the values of B-^, Bg> and B^ are plotted against cd s /o> c 

for — * i, 2, and 5 . Equation (All) may "be written simply 
m 2 



(A12) 
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Force in spring .- The force in the spring after landing is the spring 
stiffness times the displacement at the position x=0 . By use of 
equation (a 8) , the force is found to "be 


Sv 


t 2 2(l + cos 9 cosh 9 \ 

-h 1 £ ~ V i (“3) 


U=0) 


P -1^2 


e n* 
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APPENDIX B 

APPBOXIMATE SOLUTIONS 
Stresses Computed for Eigid Wing 


Method A - based on acceleration obtained with rigid wing (no 
structural" elasticity ) . - If the wings of the structure shown in figure 1 
were rigid the landing operation would he simply that of a rigid mass 
equal to. M + m alighting on a spring. The motion after arrest would he 
that of a simple oscillator having a mass M + m and a maximum velocity v. 
The solution based on these assumptions is designated method A. The 
maximum acceleration for such an oscillator is 


a = 



V CDg 


CB1) 


The bending stress is computed on the basis that the wing is loaded with 
a uniform load having an intensity equal to the mass per unit length times 
the maximum acceleration. Prom the static theory of the bending of a 
cantilever beam, the bending moment at the root resulting from this loading 
would be 


M = a — 

L 2 2 


The bending stress due to this bending moment is 


My _ amLy _ 

I 21 21 


(B2) 


With the use 

of the notation c 2 = — 

7 

frequency of 

a beam, cp = pc ( 

c I 2 

be written 



CD 


B ^ 4 


= A — — E 


— a E 


c p 


(B3) 
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This equation Is of the* same form as equation (A9) . It is noted 
that the- stress varies linearly with a> s /o5 c , and has no explicit- 
dependence on the ratio M/m. The value of A obtained by method A 
is plotted against 0) 8 /a> c in figure 3. 

Method B - based on reaction obtained with elastic wing .- In method B, 
a force equal to the maximum force given by equation (AI3) Is applied at 
the root. The rigid wing and fuselage mass then have an acceleration 


a = F -BS£- 
M + m 


With this acceleration applied to the root^ the static stresses induced 
at the root would be 


a = — 
I 


amLy 

21 

3 

v y ml S 
c p 21 M + m 


co 2 


4^ 


'0 2(1 + cos e n cosh 9 n ) 


V A 


sin o^t jf 


max 


= E - — 
c p 


w 

I 


“0%=, 2 0 n 2 ( 1 + 008 e n cosh s n) 


sin cc^t j_ 


n=l ^ 


max 


= e ns 

c p 


m 


The value of A obtained by method B is plotted against cDg/o> 0 for 


M 

for — = 2 in figure 3 , 
m 


Method C - based on acceleration obtained with elastic wing .- In 
method C, the root of the rigid wing is given an acceleration equal to 
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the maximum acceleration given' by equation (A12) . The "bending stresses 
at the root is then 


cr = ^ 
I 


21 


v y ml/ 
c p 21 


3 


OO 

I 

n=l 


2 F(0) 


®o' 7F” 8ln V 1 
9 n A 


max 


E II 
c p 


Si 


2 2 

a>o 0 ( 1 + cos 0_ cosh 0 ) 
u n\ n n/ 


A 


sin a) t i 
n 


max 


E — — A 
c p 


(B5) 


The value of A obtained by method C is plotted against m a /m c 

for — = 2 in figure 3 • 
m 


Stresses Computed for Elastic Wing 


Method D - based on reaction obtained -with rigid wing .- The reaction 
that results from landing when the wing is rigid will vary sinusoidally 


with an amplitude S — , thus 

m s 


H 



sin ojgt 



8 


If , In accordance with method D, this reaction were suddenly applied to the root of the elastic 
wing, tending vibrations would he set up in the wing. The response can he found in a ma nn er 
similar to that used in the exact solution of appendix A- The only difference is that the last 
boundary condition is changed. The last boundary condition, for this case is 

m (S) +h (§) - - 8 £ i 

x=0 x=0 


With this and the remaining boundary conditions the deflections and bending stresses are then found 
as In the case of the exact solution. The maximum bending stress for this case is found to be 



In this expression, 0 n is the nth positive root of the equation 


M / \ 

1 + cos 0 cosh 0 + — (cosh 0 sin 0 + slnh 0 n cos 9 U J = 0 

ffi0TL \ 
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If A is used to replace the "bracketed term in equation (B6) 


a 


= aHe 

c p 


The value of 5 obtained "by method D is plotted against 
for ^ = 2 in figure 6. 

Method E - "based on acceleration obtained .with rigid wing .- The 
accelerations that result from landing when the wing is rigid will vary 
sinusoidally with an amplitude given by equation (Bl) , thus 


a = v ffi g sin oigt 


If, in accordance with method E, the root of the elastic wing were 
suddenly given an acceleration character! zed by the equation, the wing 
would be set into bending oscillations. Again, the oscillations or 
response can be found in a manner similar to that used in computing the 
landing response by the exact solution given in appendix A. The last 
boundary condition for this case is, however. 



V Oig Sin CJJgt / 


With this and the. remaining boundary conditions, the deflections and 
bending stresses are then found as in the case of the exact solution. 
The maximum root bending stress is found to be 


CT = I - £ 
C p 


E 0 sin CD s t j_ + 



F n sin CD n t 1 


(B7) 


n=l 


- 1 


“n 


(max 
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where 


v = 


sin 9q sinh 0 Q 

1 + cos cosh 0„ 

0 0 


F = 


sin 0 n sinh 0 n 


n 0 n fsin 0 n cosh 0 n - cos 0 n sinh 0 n ^ 


e 0 - \l 3-52 — 

“o 


0 n = nth positive root of 1 + cos 0 n cosh 0 n = 0 


'n 




Vm 


s Vm + m 

9 2 N 

- pc 

Xj 


If A is used to replace the bracketed term in equation (B7) 


— v v 

CT = A - — E 
c p 


The value of A obtained by method E is plotted against <x> B /a> c 


for — = 2 in figure 6. 
m 
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Method F - Blot and Bisplinghoff method .- In the statistical 
approach suggested hy Biot and Bispllnghoff (reference 2) , the maximum 
force and duration of impact have to he known. This method (designated 
method F herein) was applied to the case considered hy taking the 
maximum force equal to the maximum value given hy the exact solution 
(equation (AI 3 ) ) and hy taking the vertical impulse period Tj equal to 

one-half the natural period of the airplane structure with rigid wings and 


with the bottom of the spring fixed in positionj thus Tj = ~ T g , 
where 


rp _ 2 It 

s " <z 

s 


The ratio of the impulse period to the period of the nth mode — 

^n 

is found hy the following consideration. The natural frequency of the 
free -free modes of the structure , which are the modes used in the Biot 
and Bispllnghoff method, are found from the equation for frequency given 
in the exact solution (appendix A) : 


- PC 



where 0 n is taken to correspond to a structure without landing springs 

qv3 

and may he taken from figure J at - 0. With the use of the equation 

FI 

a > 0 = pc 3* - } can may he written 
Tj 


for frequency of a cantilever. 


®n " “C3-J2 


Division through hy u> s and use of the relation between period and 
frequency result in the relations 


_ ^c 9 n _ T S 
“a “ s 3-52 T n 
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Since T g is taken as 2 Tt> the equation for — .may he written 

i n 

directly 

= 1 _ 1 e n 2 

T n “ 2 as B /o 3 0 3.52 

For the ratio — = 2, the value of 9 is 4.00 and 9 ^ is 23-00. 
m c 

The value of 9 ^ is not given because it was found that the third mode 
could he neglected. 

In computing stresses, the response factor was taken directly from 
the envelope curve given in figure 13 of reference 4 . The values of 
stress coefficient obtained by method F are shown in figure 6 . 
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Figure 3.- Bending -stress coefficients when structural elasticity is 
neglected in computing stresses (that is, in stage 2). M = 2. 



Figure 4.- Physical assumptions when structural elasticity is 
considered in determining reaction or acceleration but 
neglected in computing stresses . 
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Figure 5.- Physical assumptions when structural elasticity is 
neglected in determining reaction or acceleration but is 
considered in computing stresses. 
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Figure 6.- Bending -stress coefficients when structural elasticity is 
neglected in computing reaction or acceleration (that is, in 

stage 1). M = 2. a = A-^E. - • 
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